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AERONAUTIC SYMBOLS 

1. FUNDAMENTAL AND DERIVED UNITS 




Metric 


English 

y 



Symbol 

Unit 

Abbrevia- 

tion 

Unit 

Abbreviation 


1 

meter^ 

. . m 

foot (or mile) . . 

ft (or mi) 

Ti-mp 

t 


s 

second (or hour) 

sec (or hr) 

Force 

F 

weight of 1 kilogram 

kg 

weight of 1 pound 

lb 


P 

horsepower (metric) _ _ 


horsepower. . 

S.1. 

fps 

Speed. 

V 

/kilometers per hour. 
\meters per second 

kph 

mps 

miles per hour. 

feet per second 


2. GENERAL SYMBOLS 


Weight==»i^ 

Standard acceleration of gravity =9.80665 m/s^ 
or 32.1740 ft/sec® 

W 

Mass=— 

9 - .. . 

Moment of inertia =mP. (Indicate axis of 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 


V Kinematic viscosity 
p . Density (mass per unit volume) 

Standard density of dry air,- 0.12497 kg-m'^-s** at 15° C 
and 760 mm; or 0.002378 lb-ft“^ sec^ 


Specific weight of 
0.07651 Ib/cu ft 


‘standard” air, 1.2255 kg/m® or 


Area 

:Area of wing 
Gap 
Span 
IChord 

Aspect ratio, ^ 

True air speed 
'Dynamic pressme, 2 pV^ 

L 

Lift, absolute coefficient 

Drag, absolute coefficient Dd=^ 

Profile drag, absolute coefficient Cdo= 


iu 

it 

Q 

Q 

R 


■qS 

Induced drag, absolute coefficient 

D, 

qS 


Parasite drag, absolute coeffident Cop 


3. AERODYNAMIC SYMBOLS 

Angl e of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds number, p ^ where lm& linear dimen- 
sion (e.g., for an airfoil of 1.0 ft chord, 100 
mph, standard pressure at 15° C, the corre- 
sponding Reynolds nmnber is 935,400; or for 
an airfoil of 1.0 m chord, 100 mps, the corre- 
sponding Reynolds number is 6,865,000) 
a Angle of attack 

6 Angle of downwash 

oo Angle of attack, infin ite aspect ratio 

at Angle of attack, induced 

aa Angle of attack, absolute (measured from zero- 

lift position) 

7 Ffight-path angle 


G 


Cross-wind force, absolute coefficient 
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INTERFERENCE METHOD FOR OBTAINING THE POTENTIAL FLOW PAST AN 

ARBITRARY CASCADE OF AIRFOILS 

By S. Katzopf, Robert S. Finn, and Jambs C. Laurence 


SUMMARY 

A procedure is presented for obtaining the pressure distribu- 
tion on an arbitrary airfoil section in cascade in a two- 
dimensional, incompressible, and nonviscous flow. The method 
considers directly the influence on a given airfoil of the rest of 
the cascade and evaluates this interference by an iterative process, 
which appeared to converge rapidly in the cases tried {about unit 
solidity, stagger angles of 0° and 45°). Two variations of the 
basic interference calculations are described. One, which is 
accurate enough for most purposes, involves the substitution of 
sources, sinks, and vortices for the interfering airfoils; the other, 
which may be desirable for the final approximation, involves a 
contour integration. The computations are simplified by the 
use of a chart presented by Betz in a related paper. The numer- 
ical labor involved, while considerable, is less than that required 
by the present methods of conformal transformation. Illus- 
trative examples are included. 

INTRODUCTION 

The rapid increase of interest in the design of fans and 
turbines has led to many studies of the two-dimensional flow 
past infinite lattices. Most of these studies involve approxi- 
mate procedures (for example, references 1 to 3) or present 
solutions for special classes of shapes (references 4 and 5). 
Recently, attempts have been made to obtain exact solutions 
by conformal transformation of the lattice to a circle. To 
this end, Howell (reference 6) used a procedure that first 
transformed the lattice to an isolated S-shape figure, which 
could then be transformed to a near-circle by successive 
Joukowski transformations and finally to a chcle by the 
method of reference 7. In reference 8 the cascade was 
transformed first to a near-circle and then to a circle, also 
with the use of several stages of conformal mapping. In 
reference 9 the lattice was mapped into a lattice of straight 
parallel lines by means of a function that was determined 
with the aid of the transformation of this line lattice to a 
circle. (See references 10 and 11.) These transformations 
are of considerable interest, theoretically. The methods of 
references 6 and 8 require lengthy computations, however, 
and difficulty has been experienced in obtaining accm’ate 
numerical results with the method of reference 9. All three 
methods require modifications for highly cambered contours 
or for lattices of high stagger and solidity. 

The method presented herein does not seek a conformal 
transformation directly but, like the older approximate 
methods, seeks to evaluate the interference at each airfoil 
due to the presence of all the other airfoils of the cascade,. 


The velocity distribution on each airfoil is considered to be 
the sum of that corresponding to its presence in the uniform 
free-stream flow plus that corresponding to its presence in 
the interference flow. The interference is calculated from 
the velocity distribution on the airfoils so that the method 
reduces to an iteration process in which, for the first approxi- 
mation, the interference is computed by assuming the free- 
stream velocity distribution to exist on each airfoil, and in 
subsequent approximations this velocity is corrected accord- 
ing to the interference derived in the preceding approxi- 
mation. A solution is thus found for an arbitrarily specified 
angle of attack, and this solution is used to find the conformal 
transformation to the circle and thence the solution for any 
other angle of attack. 

The present method has been found appreciably less 
laborious than the methods that seek the conformal trans- 
formation directly and is also considered more flexible in 
that it may be adapted to a variety of cascade problems that 
would be difficult to solve by formal transformation methods; 
for example, the problem of the flow about double cascades 
(or superimposed lattices) or certain types of “inverse” 
problems involving the determination of the setting or 
solidity for a given airfoil in cascade. Some of the features 
of the interference and iteration methods used should also 
be useful in the solution of flows involving a finite number of 
interfering bodies. 

SYMBOLS 

W flow function (complex potential) 

$ velocity potential 

^ stream function 

V velocity at infinity 

r circulation 

K mapping-function parameter 

V local velocity 

V vortex strength 

m source strength 

z complex variable of physical plane {x-\-iy) 
z' fixed point in physical plane 

f complex variable of reference plane (^fl-fi?) 

c profile chord 

Cl profile chord used in transformation of referenee 7 
d cascade spacing (distance between corresponding 

points on adjacent blades; see fig. 1) 

<p central angle of perfect circle obtained in trans- 
formation of reference 7 
B central angle of unit circle of figure 1 
s surface length on profile 
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^ -plane 



L 



Figure Flow singularities in circle piano and corresponding velocity vectors 
in physical plane. 

)S blade angle (angle between stagger line and normal 
to chords; see fig. 1) 

a solidity (ratio of chord to distance between profiles) 

X angle between flow direction and normal to stagger 

line 

a angle of attack relative to blade chord 

angle of zero lift for cascade, relative to blade chord 
static pressure rise 
CO turning angle of flow 

p density of fluid 


a„, Fourier series coefficients 

Subscripts: 

/ free stream 

d disturbance 

c compensating 

r due to circulation change 

ffi additional 

T total 

t tail stagnation point 

n nose stagnation point 

TE trailing edge 

s due to source rows 

® due to vortex rows 

z physical plane 

f reference plane 

0 mean flow 

1 incoming flow 

2 outgoing flow 

Xo at flow direction Xo 

Xo' at flow direction Xo' 

THEORY OF INTEEFEEENCE CALCULATIONS 

In order to explain better the basic concepts and pro- 
cedures of the interference calculations, discussion of the 
iteration steps will be postponed for the present, and the 
interference calculations will be described as if they were 
being used to verify a known solution. 

Breakdown of the flow function into four components. — 
Attention is fixed on one airfoil of the infinite cascade which 
will be designated the central airfoil. The flow function on 
the boundary of this airfoil is considered to be the sum of 
the following components : 

Wf the flow function for the central airfoil, considered as 
isolated in the free-stream flow (the vector average 
of the flow far in front of the cascade and the flow 
far behind the cascade) . Inasmuch as the boundary 
is a streamline in this flow, Wf=^f. 

Wi the disturbance along the contour caused by the pres- 
ence of aU the other au’foils of the cascade, designated 
the external airfoils {Wa=^a+'i'^a)- 
Wc the compensating flow function (which may have 
singularities only within the central airfoil) that is 
required to maintain the airfoil a streamhne in the 
presence of the distmbance flow. It is determined 
by the condition that, on the boundary, its stream 
function must be equal and opposite to the distm’b- 
ance stream function. Thus, IFc=#e+'i^c, where 

Wr the contribution of the circulation that must be added 
to maintain the traihng-edge condition; it has only a 
real component (IFr=$r). 

The sum Wa-^Wc+Wv represents the net change of flow 
function due to the presence of the external airfoils; it wiU be 
designated the additional flow function Wa=^a- The sum 
Wa-{-Wf will be designated the total flow function Wt=^t- 

The evaluation of the isolated, or free-stream, flow is 
readily performed by the method of reference 7 and requires 
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no further discussion in the present paper. The disturbance 
flow can be calculated when the potential distribution (or 
velocity distribution) on the external airfoils is known. 
Finally, the compensating flow and the circulation flow are 
readily determined, as will be shown, when the disturbance 
flow is known. In the following sections two methods of 
calculating the disturbance flow will be described: the 
approximate source-vortex method and the exact contour- 
integral method. 

Disturbance flow by approximate source-vortex method. — 

Each of the external airfoils is considered to be adequately 
represented by an arrangement of about two sources, three 
sinks (or negative sources) , and five vortices distributed along 
its mean line. The strengths and locations of these singulari- 
ties are chosen on the basis of the chordwise thickness dis- 
tribution and chordwise velocity distribution. The choice is 
somewhat arbitrary and may be left to the judgment of the 
worker; however, a detailed method of choice has been 
described in the section entitled “Computational Methods.” 
The disturbance flow, then, is that of about ten infinite rows 
of singularities, equally spaced along the cascade direction 
except that none are located where the central ahfoil is to be 
placed. The field of each vortex row is shown in figure 2 
where, for convenience, the vortices are assumed to be of 
unit strength, spaced at unit distance along the y-axis. 
This figure is from reference 1 and the equation for the flow 
is (reference 2) 

loge sinh TTS+tT- loge TTZ 
Ztt 

In order to find the contribution to the disturbance flow 
caused by a row of vortices at, say, 0.3 chord on the external 
airfoils, the central airfoil, drawn to scale and properly 
oriented relative to the cascade direction, is placed at the 
center of figure 2, with the origin at 0.3 chord on the mean 
line. The values of velocity potential and stream function 
read at selected points along the airfoil contour, multiplied 
by the assumed vortex strength, give directly the contribu- 
tion of this vortex row to 9a and 'i'a. By shifting the central 
airfoil so that the origin is located, in turn, at each of the 
other assumed vortex positions along the mean line and 
repeating the foregoing process, the contributions of all the 
vortices in the external airfoils are obtained at the same 
points. The sum of these values at a given point on the 
central airfoil represents the contribution of the vortices in 
the lattice to the disturbance function Wa at that point. 
The contributions of the sources can be found in the same 
way except that the lines marked are considered as — $ 
and the lines marked $ are considered as Sinks are 
considered as negative sources. 

Contour-integral method for evaluating disturbance flow 
function. — In the preceding section, the disturbance field 
was calculated approximately by representing each airfoil by 
a somewhat arbitrary arrangement of vortices, sources, and 
sinks distributed on the mean line. An airfoil may be 
represented exactly by a continuous distribution of vortices 
along its contour, the linear density of which at every point 
equals the velocity on the airfoil at that point (reference 12). 
The field at a point on the central airfoil due to a row of 


corresponding surface elements of the external airfoils (that 
is, a row of vortices of strength % ds) may be obtained 
directly from figure 2. Integration of this contribution 
along the contours of the external airfoils provides an exact 
determination of the disturbance field. The procedure is an 
obvious modification of the preceding approximate method. 

Let 9 and 'S' (without subscripts) denote, respectively, the 
potential and stream function of the row of unit vortices in 
figure 2. In order to determine the disturbance potential 
and stream function at a point s' on the central airfoil, the 
airfoil contour, drawn to scale and correctly oriented relative 
to the cascade direction, is superimposed on figure 2 so that 
the origin falls, in turn, at a number of points s on the con- 
tour, and for each setting values of 9 and ^ are read at the 
point s'. Then the disturbance flow function at s' is given 

by 

4>d(s') = J^4^3 '(s) ds 

'9a(z') = ds 

where 

Vt(.z) local velocity on the airfoil at variable point s 
s distance along airfoil contour 
$, ’4' values read at s' when origin of figure 2 is at s 
and the integration is performed along the airfoil contour. 
Since Vj-(z) ds=d9r{z), the foregoing equations can be re- 
written as 

^i(2') = d9r{z) 

~9a{z') = ^9 d9r{z) 

so that the disturbance potential and stream function at 
point s' are readily evaluated by plotting 4> and 9 against 
4>j. and measuring the area under the curv^es. 

Determination of compensating flow and circulation 
flow. — As has been indicated, the compensating flow func- 
tion may have singularities only within the central airfoil 
contom", and on the contour, the stream function must be 
exactly equal and opposite to the disturbance stream func- 
tion. From the known transformation of the isolated airfoil 
to the circle, which was found in the process of determining 
Wf, the correspondence between points on the airfoil and 
points on the circle is known. If, then, the desired com- 
pensating stream function is expanded as a Fourier series 
in terms of the circle angle >p, 

CD 

(®n COS 6„ sin n<p) 

71 = 1 

its corresponding velocity potential will be (reference 7) the 
conjugate series 

(— cos 7i(p-\-a„ sin n<p) 

n=\ 

The determination of 9c from 9c is readily accomplished by 
the method of reference 13. 
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In order to maintain the trailing-edge condition, a vortex 
Fa must be added at the center of the circle of such strength 
that Val2it equals the value of —d^Jdip at the trailing edge 
(determined numerically or graphically from a faired plot 
of against <p). The corresponding contribution to the 
potential is 


The velocity potential that constitutes 

the net effect of putting the airfod in the cascade (that is, 
the net interference effect) may now be determined by simple 
addition of the three components. Presumably, since the 
calculations were made with the correct should be the 

difference between and 3>i.. 

In the final step, is differentiated with respect to dis- 
tance along the airfoil to get the corresponding interference 
effect Va on the velocity, which should be the difference be- 
tween Vf and Vj.. Convenient procedures for performing 
these calculations are discussed in the section entitled “Com- 
putational Methods.” 

ITERATION METHOD 

In the preceding sections the basic concepts and procedures 
of cascade interference calculations have been outlined. In 
the present section, the application of such calculations in 
the proposed iteration method of solving cascade flow will 
be discussed. 

As first attempted, the method was essentially as follows: 
In the first step, is assumed to equal <I>/and a first approxi- 
mation to is calculated on this basis by the methods just 
described. In the second step, is assumed equal to the 
sum of and this first approximation to and a second 
approximation to is computed. The succeeding steps 
follow the same pattern and are continued until two succes- 
sive distributions are essentially the same. The source- 
vortex method was used for the earlier approximations, but 
the final approximation, when convergence is practically 
complete, was made by the contour-integral method. This 
procedure, however, was found to converge relatively slowly 
in some cases, and the general practicability of the interfer- 
ence method depends on a slight modification of the source- 
vortex method. 

The modification depends upon the observations that the 
contribution of the sources and sinks to changes by rel- 
atively little from one approximation to the next and that 
the contribution of the vortices to <i>a is nearly proportional 
to their total strength and relatively independent of then- 
distribution. Obviously, if it were exactly true that the 
contribution of the sources and sinks is constant and that 
the contribution of the vortices is proportional to their total 
strength, only one interference calculation would be required 
and the solution could then be obtained through a simple 
algebraic equation. Thus, let 


Fa^ contribution of vortices to Fa when F/ is assumed on all 
external airfoils 

Then, by the preceding assumptions. 


whence 



( 1 ) 


Since the assumptions are not exactly true, the value of 
Fy so calculated is correspondingly inexact; however, it is 
much closer to the true value than if it were taken simply 
as Fi-l-Faj-f-Faj^. Correspondingly, the potential 


is much more accurate than the sum 

The second approximation is similarly adjusted. Thus, 
corresponding to the 4>r distribution just obtained, a new 
set of sources, sinks, and vortices are distributed along the 
mean line, and new values of F^^ and Fa„ are calculated. 
Adjustment follows, as before, from the equation 

rr,= r,-i-r.^+^r„„ 


where the subscripts 1 and 2 refer to the first and second 
approximations, respectively. Solving for Fj.^ gives 


1 

Tt. 


( 2 ) 


and, finally, the potential is given by 

This simple modification of the procedure is so effective 
that, in the cases tried, the first step gave solutions that 
would be satisfactory for many purposes and the procedure 
had practically converged at the second step. The addi- 
tional complication of keeping the source-sink and the 
vortex effects separate so that F^^ and Fa„ can be separately 
computed is relatively minor and amply repaid by the 
rapidity of convergence. 

After the source-vortex method has essentially converged, 
a final approximation by the contour-integral method may 
be desirable. In the cases computed, however, this final 
step was found to introduce only minor changes in the result. 

THE FLOW AT OTHER ANGLES OF ATTACK 


Fr total circulation on airfoil in cascade 

F/ total circulation on isolated airfoil at same angle of attack 

Fa additional circulation (Fg.— F^) 

Ta, constant contribution of sources and sinks to F„ 


From a known velocity distribution at a given angle of 
attack, the angle of zero lift and the slope of the lift curve, 
together with the velocity distribution at any other angle 
of attack, may be obtained. For this purpose, the lattice 


IL 
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is conveniently considered to be related conformally to an 
isolated circle by a periodic transformation, which might 
be, say, of the type used in reference 6, 8, or 9. The explicit 
form of the transformation, however, is not needed for the 
present purpose. 

The flow function in the circle (f) plane that corresponds 
to the desired flow in the physical ( 2 ) plane is 


W= 



f loge 



f2_g2«: 

^2_e~2K 


(3) 


In the f-plane this flow may be interpreted as that due to 
the system of sources, sinks, and vortices shown in figure 1. 
The unit circle is a streamline of the flow and the 

circulation about any contour enclosing this circle but not 
enclosing the points f— is r (positive clockwise). 

In the physical ( 2 ) plane, the complex velocities at the 
points 2 = 0 = and 2 =— are determined by equation (3) and 
the transformation. Thus, 

and 

where the angles and velocities are defined in figure 1. The 
flow far from the lattice is seen to be the same as that of an 
infinite vortex row in the uniform flow —Voe~‘“K It should 
be noted (fig. 1) that \o=ao+fi, Xi = ai + /3, and \2=o!2+i3. 
In the following paragraphs it will be shown how to obtain 
from the given solution in cascade the parameter K and the 
stagnation points and 6t for the corresponding flow about 
the circle. These values fix the angle of zero lift and the 
slope of the lift curve of the airfoil in cascade; together with 
the known potential distribution they determine the con- 
formal correspondence between the profile and the circle 
and, hence, the velocity distribution at any angle of attack. 

Since the airfoil contour ( 2 -plane) is conformally related 
to the unit circle (f-plane), it follows that at any given angle 
of attack «o the change of velocity potential from nose to 
tail stagnation point on both upper and lower surfaces must 
be the same for the circle and for the profile in cascade. 
These potential changes can readily be obtained for the 
single solution on the lattice from the final <I>r distribution. 


The velocity potential on the unit circle is obtained from 
equation (3). Thus, 


cFo 


1 

2x0- 


^cos Xo logc ^ 


cosh jjC-l-cos 6 \ 
cosh K—cos 6/ 


+ 


2 sin Xo tan ' 


sin 0 , r , tan 0 ”] 

sinhir^Fod tanhKj 


(4) 


and the change of potential from nose stagnation point to 
tail stagnation point 0i is 


cF 


| ~(cosh K—cos gJ(cosh fjC+cos 0„) 
cos 0,)(cosh if— cos 0„) 


>f If , r (cosh K— 

2 sin X, tan- r 

1_ smmit-1-sin sm 0, J 

_i p (tan tan 0,) tanh .K Jj 


] 


+ 


Fod 


tan“ 


[_ tanh^ii'-j-tan 6^ tan dt 


+ 


(5) 


This potential change may be obtained for either the upper 
or the lower surface. Two values are obtained depending 
on the choice of quadrant for the third term of equation (5). 
The condition of zero velocity at nose and tail stagnation 
points is 

p 

sin 0 cos Xq— cos 0 tanh K sin ^o~ 2 ~ y g K=0 (6) 


By use of the known values of F, A>I>r, and Xo, equations (5) 
and (6) can be solved simultaneously for 0t, and K. 
Equation (6) can be considered as a quadratic in sin 0 and 
with an assumed value of K determines corresponding values 
of 0n. and 0f Equation (5) then determines A$f. By the 
proper choice of values of K, a curve of A$f against K may 
be plotted such that at a point on this curve A$f=A$r. 
The value of K at this point is the desired value; the corre- 
sponding values of 0n and 6, are then given by equation (6). 
A convenient initial choice for K is the value that corresponds 
to a lattice of straight lines of the same stagger and of about 
10 percent or 20 percent higher solidity. Figure 3 is of aid 
in this respect. The computed values of K and together 
with equation (6), determine the angle of zero lift (F=0) 
with respect to the airfoil chord, thus 


jj = tan~^ 


tan 0t 
tanh K 




(7) 


and the slope of the lift curve, based on mean velocity, is 
obtained by differentiating equation (6) with respect toXo; 
thus, 

dc;_4 Vsin^ g,-|-sinb^ K ,o\ 

doco a sinh K cosh K 


A correspondence between points on the airfoil and points on the unit circle may be obtained by comparing the values 
of computed by equation (4) with the values of from the known potential distribution. The points {x,y) on the 
profile and 0 on the circle for which $!-=4>r are corresponding points. The velocity on the lattice profile for the stream 
angle Xo is 


^ = 


1 


cos Xo cosh K (sin sin 0()~sin Xo sinh K (cos 0— cos 0^~ 


dz 

Vo~d0~ 

dz ir 

cosh^Ai — cos^ e 
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where the term in brackets, which represents the velocity on the circle boundary, is obtained by differentiating equa- 
tion (4). It follows that the velocity corresponding to a new stream angle Xo' is 


f V \ _/v\ r cos Xq' cosh K (sin g— sin g,)— sin Xp' sinh K (cos g— cos d,) 
\T^oAo'~\T^AoLcos Xo cosh K (sin 0— sin fl,)— sin Xo sinh K (cos 0— cos dt) _ 



0 15 30 45 60 75 90 

Blade angle^ /3, deg 


Figure 3.— Relation between solidity, blade angle, and parameter K for an infinite lattice of 

flat plates. 


The following relations, which describe the flow far away 
from the lattice, are of interest. The stream angles X, and 
X 2 at and 2 = — °o are 


and 


Xi = tan”* 


sin Xod 
cos 


r 

~2Uod 

Xo 


X2=tan“* 


sin Xo 


cos 


r 

Xo 


and the angle through which the fluid is turned by the lat- 
tice is given by 



The rise in static pressure across the lattice is 



= cos" Xo (sec^ Xi — sec^ X 2 ) 

REMARKS ON CONTOUR MODIFICATIONS CORRESPONDING 
TO LOCAL PRESSURE CHANGES 

In reference 14, the modification of an airfoil contour to 
obtain, approximately, desired small changes in the pressure 
distribution is discussed. The method, based on the for- 
mulas of reference 7, evaluates a slight modification of the 
conformal transformation of the circle to the airfoil, such that 
the stretching factor at every point is changed in proportion 
to the. desired relative change in local velocity. 

Although in reference 14 the airfoil was assumed to lie in 
a straight uniform field, the treatment is equally applicable 
when the airfoil is in a curved or distorted flow field. Accord- 
ingly, the procedure should be applicable to airfoils in cascade, 
provided the same modification of the external airfoils leaves 
the disturbance flow field essentially unaffected. This 
condition may not always be satisfied; however, in such 
cases the method could possibly be improved by a procedure 
analogous to that described in the section of the present 
paper entitled “Iteration Method.” 

COMPUTATIONAL METHODS 

The basic theory has been presented. In the following 
sections some of the methods used for performing the actual 
computations will be discussed. 

Selection of points for evaluation of disturbance flow. — 
The determination of the compensating flow by the method 
of reference 13 requires that the disturbance flow be evaluated 
at points that, by the conformal transformation of the 
isolated airfoil to a circle, correspond to points equally 
spaced about the circle. These points, which are located by 
reference to the conformal transformation, are preferably 
chosen so that one is at the trailing edge. Experience has 
shown that 12 points at 30° intervals yield acceptable results, 
at least in the first approximation, but that 24 points at 15° 
intervals are desirable in the final approximation for best 
accuracy. If the final approximation is by the contour- 
integral method, which is rather tedious, an acceptable 
compromise is as follows: Evaluate and 'Pi directly at 15° 
intervals only in the region of the leading edge and at 30° 
intervals over the remainder of the contour, interpolating 
from a faired curve to determine and for the remain- 
ing 15° interval points. 
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Inasmuch as values for the 12-point and 24-point methods 
are not included in reference 13 the following table is pre- 
sented: 


k 

Cic 

7i=6 

71=12 

1 

0. 62201 

0. 63298 

3 

. 16667 

.20118 

5 

.04466 

.10860 

7 


.06394 

9 


. 03452 

11 


.01097 


Evaluation of and Va- — Integration of equation (36) of 
reference 7 along the circle boundary yields the values of the 
potential 4>/ at points on the airfoil as follows: 

^=2ae^“[ss sin (a-fd) — cos (a-|-(p)] (10) 

V 0 

where 

a angle of attack 

/3 angle of attack for zero lift 

radius of the circle to which the airfoil transforms 
<p angular position along the circle, as determined by 

the transformation 

If the transformation has been performed as recommended 
in reference 7, the constant a will be slightly less than one- 
fourth the chord. Altbough the potential discontinuity (cor- 
responding to the circulation) may, without loss of gener- 
ality, be placed at any point on the contour, the trailing 
edge will generally be found to be the most convenient 
location. 

The additional velocity Va is given by the derivative along 
the surface it may be determined by numerically or 

graphically differentiating $0 with respect to the circle 
angle <p and multiplying this slope by Thus, 

Va _ 1 1 d^a dip , . 

Vo~Vo ds Vo dip ds ^ ’ 

The value of ^ may be obtained from equations (37) and 
(38) of reference 7. Thus, 


or 




as 




2 J[l +(^)’l [sinl«+siii’«] 


1-b — 

d(p do 

2a-Y/[^ 1 + [sinh^'J' -b sin^e] 


( 12 ) 


where the symbols e, 9, and ^ are defined in reference 7. 


The cascade solidity need be taken into account only when 
the airfoil sketch to be used with figure 2 is constructed. 
For the subsequent calculations, any convenient airfoil 
chord may be used, provided only that the same chord is 
used for the external airfoils and for the central airfoil. The 
reason is as follows: The strengths of the sources, sinks, and 
vortices used to represent the external airfoils are propor- 
tional to the assumed airfoil chord; hence the additional 
potentials induced on the central airfoil will be proportional 
to the assumed chord. Since both the additional potential 
and the distance s along the contour are proportional to 

d^ 

the chord, the additional velocity will be independ- 


ent of the chord. 

The chord may then conveniently be chosen as that cor- 
responding to a = 1 since a would then not appear in equations 
(10) and (12). 

The net velocity at a point on the airfoil surface is the 
algebraic sum of the velocity on the isolated airfoil and the 
induced velocity Va at that point. 

Selection of vortices for source-vortex method. — For cas- 
cades of about unit solidity, the vortex distribution for an 
airfoil of conventional design may be represented by five 
vortices spaced on the moan line at 0.1, 0.3, 0.5, 0.7, and 
0.9 of the chord. The strengths of the vortices are deter- 
mined by the known chordwisc distribution of potential 
$7. on the upper and lower surfaces for the given approxima- 
tion. Thus, the difference in potential between tbe upper 
and lower surfaces at 0.2 chord is approximately the total 
vorticity between the leading edge and 0.2 chord and is 
considered to be concentrated in the vortex at 0.1 chord; 
similarly, the increase in this potential difference between 
0.2 chord and 0.4 chord yields the strength of the vortex at 
0.3 chord, and so on. The total vortex strength must satisfy 


the equation 


r 

eVo 2 


Selection of sources and sinks for source-vortex method. — 
The selection of sources and sinks to represent the thickness 
distribution of airfoils is less readily systematized than is the 
selection of vortices to represent the lift distribution. For 
conventional airfoils, a reasonably satisfactory representa- 
tion is generally attainable with a source at about 0.025 
chord, a second source midway between the nose and the 
position of maximum thickness, and sinks at 0.5, 0.7, and 
0.9 of the chord. The strength of each source or sink is 
taken as the difference between the “internal flow” at a 
station midway between it and the preceding source, and 
the internal flow at a station midway between it and the 
following source. This internal flow at a given station is 
estimated to be the product of the thickness and the average 
of the upper and lower surface velocities at that station. 

Obviously, not all airfoil shapes will be best treated accord- 
ing to the pattern just described; however, little ingenuity is 
required to adjust the treatment to a particular shape. 
In any case, the total source strength must equal the total 
sink strength. 
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PROCEDURE 


A suggested step-by-step procedure is as follows: 

(1) Obtain the velocities on the airfoil at the given angle 
of attack in a uniform stream by tbe method of reference 7. 
This step also determines a conformal correspondence be- 
tween points [X, y) on the airfoil and angles ^3 on a circle, 
and hence the po tential distribution 3>/by equation (10). 

(2) Using the procedure described in the section entitled 
“Computational Methods,” choose sources, sinks, and vor- 
tices to represent the airfoil. 

(3) Choose points around the airfoil at which the distm’b- 
ance function Wa is to be found; these points are conven- 
iently chosen, by reference to tbe conformal transformation, 
to correspond to 24 (or 12) equal intervals about the circle. 
By use of figure 2, determine at these points the contributions 
to "btf and of each source and vortex row. Sum separately 
the values due to sources and vortices at each point. 

(4) Form the compensating functions both for 

vortices and sources and determine the conjugate functions 

by the method of reference 13. Plot against v’ and 
determine the slope at the trailing-edge point. The relation 


r„=— 27 t 



determines the circulation changes and 


duo to the source and vortex rows. Obtain Fr l>y means 
of equation (1). 

(5) At each point 


(a) Sum the values of and due to the vortex rows 

Y 

and multiply by the ratio 

/ 

(b) Sum the values of and duo to tbe rows of 

sources and sinks. 

(c) Find $r=(r^-r^) 

(6) Sum the terms (a), (b), and (e) of step (5) to get 
plot against the circle angle (p, and determine the slopes 
at the points used in the original conformal transformation 

(step (1)) at which points the stretching factor ^ will be 

known. (Another procedure is to determine the slopes at 
the 24 equally spaced points by some numerical method and 
then to determine the stretching factors at these points by 
interpolating from the values found in the conformal trans- 
formation.) The additional velocity is given by equation 
(11); the net velocity on the airfoil surface is the sum of the 
additional velocity and the velocity on the isolated airfoil. 
The corresponding total potential is $7.=4>c+4>d+$r+4>/, 
where is known from step (1). 

Using this new potential and velocity distribution, repeat 
the procedure, starting from step (2) . The only modification 
is that Fr (step (4)) is now obtained from equation (2), and 
in step (5a) the correction factor is The process is 

continued until the changes in lift and velocity distribution 
become small. For practical purposes, the results obtained 
in this manner may be entirely satisfactory. More accurate 


results may be obtained, however, by application of the 
contour-integral method as described in the following three 
steps. 

(7) Place the airfoil drawing on figure 2 with the origin, 
in turn, at each of the 24 (or 12) points at which values are 
known from step (6) (considered as 2 -pomts), and read the 
chart at the same 24 points (considered as s'-points). As 
previously noted, some of the g'-points may be neglected. 
For each of the s'-points plot the values of # read at that 
point against the corresponding values of #r. By planimetry 
find the area between the faired emve and the #r-axis to 
determine 4>s. The value of for each point is determined 
similarly from a plot of ^ against the corresponding values 
of 


(8) Form the function determine its conjugate 

4>c; the circulation change is Fa=— 2ir and the 

\ a ^ /te 


potential 4>r = Fa 

Jtt 

(9) Sum the terms 4>c, 4>a, and >i>r to get plot against 
the circle angle <p, and measure the slopes. The velocities on 
the airfoil surface in cascade are obtained as described in 
step (6). Unless this velocity distribution differs widely 
from that obtained in the preceding approximation, it 
should not be necessary to repeat the procedure. 

The velocity distribution at another angle of attack may 
be obtained as follows; 

(a) Solve equations (5) and (6) for 0,, and K. A 
method of solution is indicated in the discussion following 
equation (6). The angle of zero lift and slope of the lift 
curve may then be obtained from equations (7) and (8). 

(b) Obtain the potential distribution $j- as a function of d 
(equation (4)); compare with the known >i>r to get a corre- 
spondence l)etween d and position on the airfoil. Equation (9) 
then yields the velocity distribution at stream angle X/. 


ILLUSTRATIVE EXAMPLES 

Example I. — The velocity distribution was obtained on the 
NACA 4412 airfoil in the configuration shown in figure 4, 
where /3=0°, (7=1.032, and Xu=9.7°. This example has boon 
treated in reference 8. In accordance with the foregoing 
procedure, results as follows were obtained: 

(1) In figure 5 is shown the chordwise velocity distribu- 
tions of the isolated airfoil at the angle of attack of 9.7°, as 
obtained in a second approximation by the method of 
reference 7. The lift coefficient at this angle of attack is 

1.67 (that is, ^^=0.837), the angle of zero lift of the airfoil 

is —4.24°, and the slope of the lift curve is 6.95 per radian. 

(2) By use of the procedure suggested in the section 
entitled “Computational Methods,” five vortices, two 
sources, and three sinks were chosen to represent the airfoil 
initially (fig. 6 and table I). 

(3) With the first location at the trailing edge, 12 locations 
on the airfoil were found corresponding to 30° intervals of 
the circle angle <p. These locations are shown in figure 6. 
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FiCfURE 4.— XACA 4412 airfoil in lattice arrangement. /S=0®; <r=1.032; Xo=9.7®. 



Figure 6.— NACA 4412 airfoil showing chosen locations of sources and vortices along mean 
line and locations at which chart readings were taken. 


(The piimed points correspond to 15° intervals.) Readings 
taken at these points from figure 2 are given in table II. 
These readings, multiplied by the appropriate source and 
vortex strengths, yielded the values of and due to 
sources and vortices given in table III. 

(4) The conjugate functions were determined by the 
12-point method and are given in table IV. The slopes of 
these functions at the trailing-edge point yielded circulation 


To Ta 

changes 0.006 and -y/=— 0.538, from which (equa- 

CVq CVo 

tion (1)) =0.513. This value corresponds to a first ap- 

proximation to the lift coefficient in cascade (c;=1.03). 

(5) In table IV are given the values of and due to 

vortex rows multiplied by the ratio - ^ (equation (1)), the 

1 / 

values of and due to source rows, and the function 


4>r=(rr~r,j 

Ztt 

(6) The additional potential 4>o=4>tf-|-‘kc+$r is plotted in 
figure 7. Slopes of this function were measured at points at 
which the stretching factor is known from step (1). The 
additional velocity v„ was then computed by equation (10); 
the algebraic sum of Va and the velocity in isolated flow 



Figure 5.— Velocities on XACA 4412 airfoil In isolated flow and in lattice arrangement. p=0®: ir=1.032; ao=9.7®. 
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O 80 120 ISO 200 240 280 320 360 

Circle anglCy deg 

Figure 7.— The additional flow function against circle angle for NACA 4412 airfoil in 
lattice arrangement. ^=0®; o-=l.032; ao=9.7®. 

yielded the cascade velocity (fig. 5). This velocity distribu- 
tion, together with the total potential formed the basis 
for a second approximation (figs. 5 and 7). Results of this 

r r 

approximation are ^^ = 0.006, ^^=—0.365, and C;=0.99. 
eVo ’ cKo 

Comparison of the velocity distribution with that of the 
first approximation shows that the process has satisfactorily 
converged. 

(7) The same 12 points around the airfoil were chosen as 

2 -points; these, together with four others at 15° intervals 
around the nose (primed points in fig. 6 ) were used as 
2 '-points. Readings from the chart (fig. 2 ) are given in 
table V. These values were plotted against total potential 
4>r (arbitrarily fixed at 0 on the lower surface at the trading 
edge). (A sample curve is shown in fig. 8 .) These curves 
were integrated by planimetry. The results — the disturb- 
ance potentials and stream functions and — are given 

in table VI. 

( 8 ) The function (table VI) was obtained by 24-point 
harmonic analysis and synthesis, with the use of interpolated 
values of for the points at which it was not found ex- 
plicitly. The slope of the curve at the trading-edge point 



Figure 8.— Typical curves for determination of and ■4's by contour-integral method. 
These curves are for point g on NACA 4412 airfoil in lattice arrangement. <r=l.Q32; 

ao=9.7°. 


yielded ^=- 


-0.344, from which a lift coefficient Cj=0.99 


was obtained. 

(9) The additional potential $a=$a-l-4>c+^r is plotted in 
figure 7. The velocity distribution was obtained as before 
and is plotted in figure 5. The process appears to have 
essentially converged. 

Simultaneous solution of equations (5) and ( 6 ) (table VII) 

to find the value of K at which gave A=: 0.3083, 

0 „= — 7.57°, and 0 , = 181. 72°. Equations (7) and ( 8 ) then 
yielded the angle of zero lift 17 =— 5.75° and the slope of the 


lift curve -v—*= 3.71. 
Otto 

ij=-5.94° and 4-= 
ttofo 


These values may be compared with 
3.71 from reference 8 . 


In figure 9 is shown a plot of the potential f>j- against d com- 
puted by equation (4) . A constant has been added to make 
the potential equal to zero on the lower surface at the trailing 
edge. The known total potential in cascade and the 
corresponding values of xjc are given in table VIII. Values 
of d, picked off the plot at points where 4>j- is equal to the 
given values of are shown in the adjacent column. The 
correspondence between airfoil position and the angle 0 is 
thus determined. For the flow angles Xo' = 1.81° and 
Xo' = — 5.94°, the velocity distributions were computed by 
equation (9) . In figure 10 these results are compared with the 
distributions given in reference 8 . The main results of the 
calculations are summarized in table IX. 

Example II. — In an effort to obtain in the simplest possible 
manner a reference solution at large blade angle, concerning 
the accuracy of which there could be little doubt, a lattice 
was derived by a modified Joukowski transformation. This 
transformation is discussed in detail in the appendix. The 
cascade configuration is shown in figure 11 where ^=45°, 
<7=1.006, and Xo=49°. This lattice will be referred to as 
the “derived airfoil lattice.” 



Figure 9. — Velocity potential on unit circle in f-plane, for NACA 4412 airfoil in lattice 
arrangement. ff-*1.032; oto=9.7®. 
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The procedure followed for the source-vortex method was 
similar to that of the first example; the calculations are 
outlined in figures 12 to 14. Because of the unusual shape 
of this profile, only one source was used and an additional 
sink was inserted at 0.3 chord (fig. 12). From a lift coeffi- 
cient Cj=0.84 in isolated flow, a single approximation yielded 
a lift coefficient C(=0.54 in cascade, which was the same as 
that derived from the solution by conformal transformation. 
Since the computed changes in vortex distribution were small, 
no further approximations were made by this method. By 
reference to the velocity distribution of this approximation 
(fig. 13), the process may be seen to have essentially con- 
verged to the correct solution. 

The final contour integration resulted in a lift coefficient 
c,=0.54 and the velocity distribution shown in figure 13. 
The main results of the calculations are summarized in 
table X. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., January 10, 1947. 



b 

Figure 12, — Derived airfoil showing chosen locations of sources and vortices along mean line 
and locations at which chart readings were taken. 
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Figure 14. — The additional flow function against circle angle for the derived airfoil lattice. 

^=45®; <7=1.006; Xo=4«®. 


I 





APPENDIX 

DERIVED AIRFOIL LATTICE 


The symbols used in the appendix are defined in figure 15 
and should not be confused with similar symbols used in the 
main text of the paper. 

Consider the transformation (reference 10), 

-i ‘os- ‘OS- 

The unit circle (f-plane) becomes a lattice of horizontal 
straight lines in the 2-plane, spaced at unit intervals along the 

7T 

stagger line, making an angle -—/3 with the axis of reals. 
The solidity of this lattice is 

2 / -1 Vsinh^jBT-f cos^^-j-cos )3 , 

v=- (cos i3 log. + 

sin /3 tan“'-pF=======') 

V smh^i?-f cos^/3/ 

This relation is plotted in figure 3. 

A closed curve enclosing the points f= ±6“*^ but not en- 
closing the points will transform by equation (Al) 

into an infinite lattice of closed shapes in the 2-plane, spaced 
in the same manner as the straight-line lattice. Such a 
curve is the circle 


— J- 

= 1.07e'«-b0.09e'3^5 


This circle, where /3=45° and fiT=0.331, becomes the lattice 
of profiles that has been referred to as the derived airfoil 
lattice. A flow for which this circle is a streamline and which, 
in the 2-plane, has no singularities outside the profiles, is that 
due to the system of sources, sinks, and vortices shown in 
figure 15. The velocity on the circle boundary due to this 
system is 



=A cos Xo+B sin Xq+C 


JL 

2Vo 


where 


, r sir 


A 

B=e~'!'o 

C=e-' 


sin (<^— 5i) 


sin («^ — 62) 


and 




H 1 — CO& (</>— 5]) 

H 2 —COS {<j>- 

”^2), 


Ji 


_i?l — cos (^ — 61) 

i?2 — cos (<^)- 

-S2), 

J, 

|_ . ...'^2 .. 
r TT / , 



5i = tan-i ^ 


r sin 5 


e^—r cos 5 


«=K”*+s) 

ml—e~’^°^|r^-\-e^‘^—2re’^ cos 6 


n 

I 



Figure 15.— Flow singularities in f-plane for derived airfoil lattice. 


m1=e~'>'<‘^|A-\-e^^A2re^ cos 5 


The constant T, which is the circulation about each profile 
(positive clockwise), is determined by the trailing-edge 
condition as 





COS Xo+^SlIl A| 


o) 


(A2) 


where A, B, and C are evaluated at the angle ^ which 
corresponds to the trailing edge of the profile. The angle 
of zero lift tj with respect to the airfoil chord is obtained from 
equation (A2) by setting T— 0; thus, 

57=— tan“*^ — 


The stretching factor from the circle to the lattice is 
dz~ 2 V -E 
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where 

Z?=(cosh 2K — cosh 2^ cos 20)^+(sinh 2^ sin 20Y 

E=4c cos^jS cosh^^L (cosh^— cos^0) + 

4 sin*®^ sinh^^ic (cosh®®-^— sin“0) — 
sin 2/3 sin 20 sinh 2K 

and i// and 0 are obtained from ^o, 4>, t, and 5 as 

. , e^o sin 4>-\-r sin 5 

e*o cos 0+r cos 5 

e*=r cos {d—S)-\--yje'^*°—r‘‘ sin^ (0—3) 

The velocity at any point on the surface of a profile is 

(Fo)r(^)Jil 
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TABLE I.— STRENGTHS OF SOURCES AND VORTICES CHOSEN TO REPRESENT THE NACA 4412 AIRFOIL LATTICE 


Vortex 
location 
(fig. 6) 

y 

cVa 

Source 

location 

(fig. 6) 

m 

cVo 

First ap- 
proxima- 
tion 

Second ap- 
proxima- 
tion 

First ap- 
proxima- 
tion 

Second ap- 
proxima- 
tion 

p 

0. 379 

0.290 

a 

0. 097 

0. 101 

d 

.184 

.098 


.044 

.044 

e 

. 128 

.062 

c 

— , 043 

— .041 

f 

.097 

.042 

1* 

045 

047 

V 

.052 

.023 

n 

-.053 

-.057 


TABLE II.— CHART READINGS FOR NACA 4412 AIRFOIL LATTICE, SOURCE- VORTEX METHOD 




4> for vortex row of unit strength 




'ir for vortex row of unit strength 


Origin at 















\ (fig. 6) 

Reading at 

(fig. 6 ) 

a 


y 

h 

« 

r 

V 

a 


y 

5 

c 

r 

*7 

a 

0. 002 

0. 008 

0. 008 

0 . oil 

0.010 

0.004 

0.001 

- 0 . 211 

- 0. 184 

- 0. 176 

- 0.120 

- 0. 070 

- 0 . 026 

- 0. 003 

b - - 

.002 

.008 

.008 

.011 

.010 

.003 

0 

—.182 

—.168 

—.150 

-.096 

—.047 

-.013 

0 

c 

.004 

.008 

.008 

.009 

.005 

-.001 

-.002 

-.116 

-.092 

-.083 

-.044 

-.011 

0 

-.010 

d 

.006 

.006 

.006 

.004 

-.002 

-.007 

-.007 

-.045 

-.030 

-.026 

-.004 

0 

-.019 

-.056 

e 

.002 

.002 

.002 

-.004 

-.012 

— . 013 

013 

-.007 

-.001 

0 

-.002 

-. 025 

-.040 

-.125 

f - 

0 

-.001 

-.002 

-.008 

-.016 

— . 013 

-. Oil 

0 

-.001 

-.001 

-.019 

-.067 

-.115 

—.174 

g - 

0 

0 

0 

-.003 

-.005 

-.002 

0 

0 

-.003 

-.002 

-.024 

-.063 

-.120 

-.183 

h 

-.002 

0 

.001 

.002 

.005 

.010 

.016 

0 

.001 

0 

-.012 

-.044 

-.090 

-.166 


-.012 

-.007 

-.006 

.001 

.008 

.015 

.023 

-, 014 

-.006 

-.004 

.001 

-.013 

-.048 

-.098 

j - 

-.020 

— . 015 

-.014 

-.006 

0 

.010 

.016 

-.058 

— . 042 

-.037 

-. 010 

.001 

010 

— . 042 

k 

-.016 

-.011 

-.010 

-.005 

.002 

0 

.004 

-.125 

-.102 

-.094 

-.060 

-.018 

0 

-.006 

1 

-.005 

0 

0 

.005 

,004 

.001 

0 

-.185 

-.161 

-.153 

-.097 

-.048 

-.015 

0 
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TABLE III.— CONTRIBUTIONS OF INDIVIDUAL SOURCE AND VORTEX ROWS TO THE DISTURBANCE FLOW FUNCTION 
ON THE NACA 4412 AIRFOIL IN CASCADE; FIRST APPROXIMATION, SOURCE-VORTEX METHOD 


Point on airfoil (flg. 6) 

Vortex location 

Source location 


5 

€ 


V 

a 

y 

c 

r 

V 

due to vortex rows 

due to source rows 


0. 0030 

0. 0020 

0. 0013 

0.0004 

0 

0. 0205 

0. 0077 

-0. 0030 

-0. 0011 

-0. 0001 

b 

.0030 

.0020 

.0013 

.0003 

0 

.0176 

.0066 

-.0020 

-.0006 

0 


.0030 

.0016 

.0006 

-.0001 

-.0001 

.0112 

.0036 

-.0005 

0 

-.0005 

d . - - 

.0023 

,0007 

-. 0002 

-.0007 

-. 0004 

.0044 

.0011 

0 

-.0008 

-.0030 


.0008 

-.0007 

-.0015 

-.0013 

-. 0007 

.0007 

0 

-.0011 

-.0018 

-.0065 

f 

-.0004 

-.0015 

-.0019 

-.0013 

-.0006 

0 

.0001 

-.0024 

-.0052 

-.0092 


0 

-.0006 

-.0006 

-.0002 

0 

0 

.0001 

-. 0027 

-.0054 

-.0097 

h - 

0 

.0005 

.0006 

.0010 

.0008 

0 

0 

-.0019 

-.0040 

-.0083 


-.0026 

.0002 

.0010 

.0014 

.0012 

.0014 

.0002 

-.0006 

-.0022 

-.0052 


-.0057 

-.0011 

0 

.0010 

.0008 

.0056 

.0016 

0 

-.0004 

-.0022 

k 

-.0042 

-.0009 

.0003 

0 

.0002 

.0121 

.0041 

-.0008 

0 

-.0003 

1 - — 

0 

.0009 

.0005 

.0001 

0 

.0180 

.0067 

-.0021 

-.0007 

0 


Si'd due to vortex rows 

Sf'rf due to source row.s 


-0. 0698 

-0. 0221 

-0.0090 

-0. 0024 

-0. 0001 

0. 0002 

0. 0004 

-0. 0004 

-0. 0002 

0 

b 

-.0599 

-.0176 

-.0060 

-.0013 

0 

.0002 

.0004 

-.0004 

-.0001 

0 


-.0348 

-.0081 

-.0014 

0 

-.0005 

.0004 

.0004 

-.0002 

0 

.0001 

d 

-.0114 

-.0007 

0 

-. 0018 

-.0029 

.0005 

.0003 

,0001 

.0003 

.0004 


-.0005 

-.0005 

-.0032 

-. 0039 

-.0065 

.0002 

.0001 

.0005 

.0006 

.0007 

f 

-.0003 

-.0035 

-.0073 

-.0111 

-.0091 

0 

.0001 

.0006 

.0006 

.0006 


-.0010 

-.0044 

-.0081 

-.0116 

-. 0096 

0 

0 

.0002 

.0001 

0 

h 

.0002 

—.0021 

-.0056 

— . 0087 

-.0081 

-.0002 

0 

-.0002 

-.0004 

-.0008 


-.0023 

-.0002 

-.0017 

-.0046 

-. 0051 

-.0012 

. -.0003 

-.0003 

-.0007 

-.0012 


-. 1057 

-.0018 

.0001 

-.0009 

-.0022 

-.0019 

-.0006 

0 

-.0004 

-.0008 

k 

-.0386 

-.0092 

-.0022 

0 

0003 

-.0016 

-.0004 

-.0001 

0 

-.0002 

1 

-.0610 

-. 0178 

-.0061 

-.0014 

0 

-.0005 

0 

-.0002 


0 


TABLE IV.— TOTAL EFFECT OF SOURCE AND VORTEX ROWS, AND CORRESPONDING DERIVED POTENTIALS AND 
VELOCITIES, ON NACA 4412 AIRFOIL IN CASCADE; FIRST APPROXIMATION, SOURCE-VORTEX METHOD 


Point 

Values at points on figure 6 

Values at points where ^ is known 

cVo 

cVq 

cVo 

cVb 

cVl 

cVo 

(-Sd+i-c) ’’’■i 
cv, r, 

4>r 

cVb 

4-0 

cVo 

cVo 


£. 

c 

1 d4*„ 
c Vo dtfi 

dtp 

^ds 

r. 

Vo 

(■) 

Vo 

O 

VT 

Vo 

o 


Sources 

Vortices 


Upper surface 


0.0000 

0. 0014 

0. 0240 

-0. 1034 

0. 0013 

0. 0067 

0. 0049 

-0. 0270 

0. 0033 

0. 6056 

0. 6089 

0. 0125- 

-0. 0356 

7.153 

-0. 254 

2. 287 

2. 033 

b 

.0001 

.0014 

.0216 

-. 0848 

.0414 

.0066 

.0295 

-.0540 

-.0015 

.5355 

.5340 

.0500 

0427 

4. 541 

-.194 

2. 002 

1.808 


.0007 

.0015 

.0138 

-.0448 

.0513 

.0050 

.0347 

—.0810 

-. 0310 

.3592 

.3282 

. 1000 

-.0456 

3. 500 

-, 160 

1. 853 

1. 693 

d 

.0016 

.0009 

.0017 

-.0168 

.0314 

.0017 

.0206 

1080 

-.0848 

.1429 

. 0581 

.2000 

-.0542 

2. 733 

-.148 

1.719 

1.571 


.0021 

-.0003 

-.0087 

-. 0146 

.0040 

-.0034 

.0004 

-. 1350 

-, 1436 

-. 0365 

-. 1801 

.4000 

-.0606 

2. 218 

-.140 

1. 523 

1.383 

f 

.0019 

-.0020 

-. 0167 

-.0313 

-.0050 

-.0057 

-.0065 

— . 1620 

1872 

-.1126 

-.2998 

.6000 

-.0506 

2. 133 

-. 108 

1.345 

1. 237 


.0003 

-.0030 

-.0177 

-.0347 

.0037 

-.0014 

.0014 

-. 1891 

— . 2084 

-.0459 

— . 2543 

.8000 

-.0226 

2. 498 

-.057 

1. 178 

1. 121 

h 

-.0016 

-.0032 

-.0142 

-.0243 

.0044 

.0029 

.0045 

-.2161 

-.2290 

.1642 

-.0648 

.9000 

-.0078 

3. 226 

-.025 

1.078 

1.053 

i 

-.0037 

-.0018 

-.0064 

-.0139 

-.0098 

.0012 

-.0053 

-.24.31 

-.2566 

.4802 

.2236 








-.0037 

.0007 

.0046 

-.0205 

-. 0356 

-.0030 

-.0236 

2701 

-.2884 

.8365 

.5481 







k 

-.0023 

.0022 

.0151 

-.0503 

-.0505 

-.0046 

-.0338 

-.2971 

31.36 

1. 1558 

.8422 



Lower surface 








— 0379! 

. om.5 

— . 0218 

—.3241 

-.3218 

1.3716 

1.0498 








!oooo 

!ooi4 

.0240 

1034 

.0013 

.0067 

.0049 

-.3511 

-. 3208 

1. 4451 

1. 1243 



















0. 0125 

-0. 0451 

8. 263 

-0. 372 

0.458 

0. 086 













.0500 

-.0651 

4. 603 

-.300 

-.314 

-. 614 













.1000 

-.0847 

3. 304 

-.280 

-. 541 

-.821 













.2000 

-.1027 

2. 413 

-. 248 

-. 676 

-.924 













.4000 

-.1157 

1.914 

-.221 

-.754 

-.975 













.6000 

-.1008 

1.903 

-.192 

-.796 

-.988 













.8000 

-.0621 

2. 335 

-.145 

-.834 

-.979 













.9000 

-. 0328 

3. 080 

-.101 

-.845 

-.946 


> Velocities along the surface are considered positive when directed from the trailing edge to the leading edge on the lower surface, and from the leading edge to the trailing edge on the upper 
surface. 
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TABLE V.— CHART READINGS FOR INTEGRATION WITH RESPECT TO NACA 4412 AIRFOIL IN CASCADE; 

CONTOUR-INTEGRAL METHOD 


■ 


a 

b 

c 

d 

e 


g 

h 

i 

i 

k 

1 


1 1.0818 

0. 6400 

0. 3106 

0.0401 

-0.2022 

-0.3280 

-0. 2911 

-0. 1064 

0.1783 

0.4986 

0. 7954 

l.OOGO 


3 0. 5876 
















0 

0 

-0. 001 

-0. 0040 

-0.0095 

— 0. 008 

0.004 

0. 020 

0. 029 

0. 0216 

0.008 

0.001 

b 

0 

0 

0 

-.003 

-.0085 

—.006 

.005 

.021 

.0275 

.020 

.0065 

0 


—.0005 

0 

0 

—.0015 

—.006 

—.005 

.006 

.0195 

.0215 

.0106 

-.0005 

-.003 

d - 

—.0045 

-.003 

-.0015 

0 

—.0016 

—.001 

.0065 

.0145 

.010 

—.0036 

-. 0125 

—.009 


—.010 

-.008 

-.006 

-.0015 

0 

.0005 

.004 

.0055 

-.006 

—.0206 

-.023 

-.016 


—.011 

—.010 

-.007 

-.002 

0 

0 

.0025 

0 

—.012 

—.025 

-.026 

-.018 

f 

—.009 

—.008 

-.005 

—.001 

.0005 

0 

.0005 

-.0025 

-.014 

-.026 

-.026 

-.014 

f' 

—.002 

—.001 

0 

.002 

.002 

0 

0 

—.0035 

—.014 

—.024 

-.0216 

-.010 


.007 

.007 

.0065 

.0065 

.0045 

0 

0 

-.002 

-.010 

-.018 

-.015 

-.002 

i' 

.015 

.015 

.012 

.011 

.0055 

0 

0 

-.0005 

-.065 

-.0125 

-.008 

.006 

h 

.024 

.022 

.019 

.0145 

.0050 

-.0035 

-.003 

0 

-.0025 

-.006 

.0005 

.015 

h' 

.028 

.027 

.022 

.014 

.001 

—.008 

—.0065 

.001 

—.0005 

—.0015 

.0055 

.0195 


.0295 

.0285 

.022 

.0095 

-.006 

—.0165 

-.0115 

.003 

0 

.0005 

.008 

.0215 

1 

.022 

.020 

.0105 

.0050 

-.0205 

-.027 

-.0185 

.006 

0 

0 

.0035 

.016 

k 

.008 

.0065 

-. 0015 

-. 012 

-.023 

-.025 

— . 014 

—.0005 

.008 

.0035 

0 

.004 


.001 

0 

-.003 

-.0085 

-.016 

-.017 

-.001 

—.015 

.021 

.0155 

.004 

0 






0 

-0. 002 

-0. 024 

-.083 

-0. 154 

-0. 205 

-0. 217 

-0. 186 

-0. 130 

-0. 066 

-0.018 

-0. 002 

b 

-.002 

0 

-.012 

-.063 

-.132 

-.178 

-.187 

—.159 

-.102 

-.045 

-.0075 

0 


-.023 

-.012 

0 

-.0205 

-.073 

-.111 

-.120 

-.095 

-.048 

-.0085 

.001 

-.0125 

d 

-.081 

-.060 

-.0205 

0 

-.018 

-.042 

-.048 

-.030 

-.004 

.0025 

-.023 

-.062 


-.153 

-. 128 

-.070 

-.0175 

0 

-.006 

-.0075 

0 

.002 

-.023 

-.077 

-.131 


-.184 

-.157 

-.091 

-.031 

-.003 

—.0005 

—.0015 

.002 

-.0035 

-.040 

100 

-.159 

f 

-.204 

-.180 

-.109 

-.042 

—.007 

0 

0 

.0015 

-.011 

-.053 

-.119 

-.181 

f 

-.214 

-.189 

-.119 

-.048 

-.0085 

0 

0 

-.001 

-.017 

-.064 

-.129 

-.191 


-.212 

-.186 

-.116 

-.047 

-.008 

0 

0 

-.0015 

-.0175 

-.063 

-.128 

-.189 


-.205 

-.176 

-.110 

-.0415 

-.005 

.0015 

0 

-.001 

-.0145 

-.056 

-.120 

-.181 

h 

-.184 

-. 158 

-.093 

-.030 

0 

.0015 

-.0015 

0 

-.008 

-.044 

-.103 

-.160 

h' 

-.158 

-. 133 

-.073 

-.017 

.0035 

-.002 

-.006 

-.002 

-.003 

-.0295 

-.082 

-.136 


-.127 

-.101 

-.048 

-.0035 

.0025 

-.011 

-.017 

-.0085 

0 

-.015 

-.058 

-.105 


-.062 

-. 044 

-.008 

.003 

-.0225 

— . 054 

-. 063 

-.044 

-.0155 

0 

-.016 

-.047 

k 

-.018 

—.0075 

.001 

-.012 

-.081 

—.123 

-.130 

-. 103 

-.059 

-.016 

0 

-.0085 

1 

-.002 

0 

-.013 

-.009 

-.131 

-.079 

-.188 

-.164 

-.108 

-.0475 

-.009 

0 






1 Upper surface at trailing edge. 

2 Lower surface at trailing edge. 


TABLE VI— DERIVED POTENTIALS AND VELOCITIES ON NACA 4412 AIRFOIL IN CASCADE; CONTOUR-INTEGRAL 

METHOD 


Values at points on figure 6 

Values where ^ is known 
ds 





4>r 

4*<i 



X 

1 

dtp 

Va 

Vf 

VT 

z' 

cVo 

cVo 

cVo 

cVq 

cVo 

cVo 

cVo 

c 

c Vo dtp 

^da 

Vo 

Fo 

V, 












Q) 

(') 

(') 

a 

0. 0281 

-0. 0687 

0. 0021 

-0.0144 

0. 0158 

0. 5717 

0. 5875 



Upper surface 



h 

.0259 

—.0571 

.0271 

— . 0432 

.0098 

.5302 

.5400 







C 

.0164 

-.0314 

.0382 

-.0720 

-.0174 

.3280 

.3106 







d 

.0033 

-.0098 

.0257 

-.1008 

-.0718 

.1119 

.0401 

0.0125 

-.0403 

7. 153 

-0. 288 

2. 287 

1.999 


—.0129 

-.0049 

.0090 

-.1296 

— . 1335 

—.0687 

—.2022 

.0500 

-.0496 

4. 541 

-.225 

2. 002 

1.777 


-.0184 

-.0071 

.0017 

-.1440 

-.1607 

-.1237 

—.2844 

. 1000 

-. 0525 

3.500 

-.184 

1.853 

1.669 

f 

-.0219 

-.0114 

-.0012 

-. 1684 

-.1815 

-.1466 

-.3280 

.2000 

-.0667 

2.733 

-.182 

1.719 

1.537 

f' 

-.0210 

-.0130 

-.0031 

1728 

-.1969 

-.1329 

-.3298 

.4000 

-.0650 

2. 218 

-.144 

1.623 

1.379 


-.0186 

-.0164 

-.0035 

-.1872 

-.2093 

-.0818 

-.2911 

.6000 

-.0452 

2. 133 

-.096 

1.346 

1.249 


-. 0157 

0151 

-.0015 

-. 2017 

-. 2189 

.0060 

-. 2129 

.8000 

-. 0189 

2. 498 

-. 047 

1.178 

1.131 

h 

-.0128 

-.0128 

-.0043 

-.2161 

-.2332 

.1268 

-.1064 

.9000 

-.0042 

3.226 

-.014 

1.078 

1.064 

W 

-.0096 

-.0118 

-.0066 

-.2305 

-.2467 

.2747 

.0280 







i 

-.0071 

-.0092 

—.0122 

-.2449 

—.2642 

.4426 

.1783 








.0002 

-.0187 

-.0261 

2737 

-.2996 

.7982 

.4986 



Lower surface 



Ir 

.0115 

—.0362 

—.0323 

— . 3025 

-.3233 

1. 1187 

.7954 







L- 

.0232 

-.0586 

-.0224 

-.3313 

-.3305 

1.3365 

1.0060 








.0281 

-.0687 

.0021 

-.3601 

-.3299 

J.4117 

1.0818 

0. 0125 

-0. 0598 

8.263 

-0. 494 

0. 458 

-0. 036 









.0500 

-.0792 

4. 603 

-.364 

-.314 

-.678 









.1000 

-.0938 

3. 304 

-.310 

-.641 

-.861 









.2000 

-.1078 

2. 413 

-.260 

-.676 

-.936 









.4000 

-.1199 

1.914 

-.230 

-.764 

-.984 









.6000 

-. 1073 

1.903 

-.203 

-.796 

-.999 









.8000 

-.0664 

2. 335 

-.132 

-.834 

-.966 









.9000 

-.0333 

3. 080 

-.103 

-.846 

-.948 


^ Velocities along the surface are considered positive when directed from the trailing edge to the leading edge on the lower surface, and from the leading edge to the trailing edge on the 
upper surface. 
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TABLE VII.— COMPUTATION OF MAPPING FUNCTION CONSTANTS FOR NACA 4412 AIRFOIL LATTICE 


K 

Bn 

(<ieg) 

Bt 

(deg) 

s 

T 

17 

cVo 

0. 3 - - 

-7.37 

181.62 

1.1307 

0. 0157 

0. 2879 

1. 4343 

0. 4 

-9. 81 

182.38 

.9585 

.0152 

.2895 

1.2632 

0. 32 

-7. 85 

181.83 

1.0920 

.0156 

.2886 

1.3962 

0.311 - - 

-7.63 

181.73 

1. 1088* 

.0157 

.2882 

1. 4127 

0. 308-- 

-7.57 

181. 72 

1.1147 

.0157 

.2882 

1.4186 

0. 309 

-7.59 

181. 73 

1.1128 

.0157 

.2882 

1.4167 

0. 3083 

-7. 57 

181. 72 

1.1144 

.0157 

.2882 



^_ CQS Xo , r (cosh K—eos 

2ira L(cosh£T+cos 


gQ(cosh JiT+cos I 


tfi)(cosh K—cos 


w-i 

wJ 


T= 


sin Xo 


U- 


r 

27tcVo 


^ , r (sin sin et) sinh K T 

L siiili2K’+sm sin dt J 

. _i r (tan gn~tan $i) tanb K "\ 
L tanhs.K'+tan Bn tan Bt J 


TABLE VIII.— RELATION BETWEEN CIRCLE ANGLE 6 (f-PLANE) AND LOCATION ON NACA 4412 AIRFOIL IN CASCADE 


Upper surface 

Lower surface 

Point 

Vo 

xfc 

e 

(deg) 

Point 

$r/c Vo 

xjc 

B 

(deg) 


-0. 2911 

0. 006 

1.4 


0. 5875 

1.000 

-178.3 

-167.3 

ii 

-. 1064 

.090 

17.0 

b 

. 5400 

.920 


.1783 

.270 

45.7 


.3106 

. 707 

—135. 8 


.4986 

.501 

98.6 

d 

.0401 

.435 

—82. 9 

k : 

.7954 

.740 

145.2 


-.2022 

. 187 

-35. 1 

1 . 

1. 0060 

.927 

168.2 

f 

-.3280 

.033 

-11.8 




TABLE IX.— CONSTANTS OF NACA 4412 AIRFOIL LATTICE 


Method 

Ar./cFo 

Ar,/cVo 

Ar«/cVo 

I Cl 

K 

Bt 

(deg) 

dci/dao 

(deg) 


0. 006 

— 0. 538 

-0.324 

1.03 






.006 

-.365 

-.342 

.99 








-.346 

.99 

0. 3083 

181. 73 

3.71 

—6. 75 

Method of reference 8 




1.00 

. 3109 1 

181. 79 J 

3.71 1 

-5.94 


TABLE X.— CONSTANTS OF DERIVED AIRFOIL LATTICE 


Method 

Ar./cVo 

AF^cTTi 

Ar«/cVo 

Cl 

K 

Bt 

(deg) 

dcijdao 

n 

(deg) 


—0.083 

—0. 101 

—0, 148 

0. 54 








—.152 

, 54 

0. 2637 

193. 50 

5. 11 

—2. 03 

Conformal transformation 




.54 

.2635 

193. 46 

5.11 

-2.11 


U. S. GOVERNMENT PRINTING OFFICE: 1949 











Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Designation 

Sym- 

bol 

Longitudinal 

Lateral 

X 

Y 

Normal 

Z 




Moment about axis 


Velocities 


Force ^ ^ " 

^«allel Linear 

Sym- Positive Designa- Sym- (compo- 
symbol Designation / . alonir Angular 


bol nent along 
axis) 


X Rolling L 

Y Pitching M 

Z Yawing. N 


Absolute coefficients of moment 
(roUmg) (pitchmg) 


(yawu^) 


Angle of set of control surface (relative to neutral 
position), S. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 

P Pow 


Diameter p p 

Geometric pitch 

Pitch ratio n a 

Inflow velocity 

Slipstream velocity a Effici 

fjv * 

Thrust, absolute coefficient Ct= — ^ Revo 

O # Effec 

Torque, absolute coefficient (7q=— 

5. NUMERICAL RELATIONS 


Power, absolute coefficient Cp— 

6 / ps 

Speed-power coefficient == ^ 
Efficiency 

Revolutions per second, rps 
Effective helix angle =tan~'^ 2 ^^^ 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 




